Suspensions Thermal Noise in the LIGO Gravitational Wave Detector 
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We present a calculation of the maximum sensitivity achievable by the LIGO Gravita- 
tional wave detector in construction, due to limiting thermal noise of its suspensions. We 
present a method to calculate thermal noise that allows the prediction of the suspension 
thermal noise in all its 6 degrees of freedom, from the energy dissipation due to the elastic- 
ity of the suspension wires. We show how this approach encompasses and explains previous 
ways to approximate the thermal noise limit in gravitational waver detectors. We show how 
this approach can be extended to more complicated suspensions to be used in future LIGO 
detectors. 
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I. INTRODUCTION 



Thermal noise is a fundamental limit to the sen- 
sitivity of gravitational wave detectors, such as the 
ones being built in the use by the LIGO project [|. 
Thermal noise is associated with sources of energy 
dissipation Q , following the Fluctuation-Dissipation 
Theorem. Thermal noise comes in at least two im- 
portant kinds: one due to the brownian motion of 
the mirrors, associated with the losses in the mir- 
rors' material; and another due to the suspension 
of the mirrors, due to the losses in the wires' ma- 
terial. The limits following from these assumptions 
(losses due to elastic properties of materials) are a 
lower limit to the noise in the detector, since there 
may always be other sources of energy dissipation in 
imperfect clamps, mirror attachments, etc. But the 
correct calculation of the thermal noise limit is es- 
sential to the design of detectors and diagnostics of 
the already-built detectors. We will deal in this arti- 
cle with thermal noise of suspensions (not of internal 
modes of the mirrors themselves), and assume only 
losses due to the elasticity of the suspension wires. 

The calculation of thermal noise can be done in 
several ways @, [|, @, §, §. All of these follow 
the Fluctuation-Dissipation Theorem (FDT), but a 
complication arises because in suspensions there are 
two sources of energy (gravitational and elastic), but 
only one of them is "lossy" (elastic energy). More- 
over, the losses in the suspension wires are associ- 
ated with their bending, and seems to be localized 
at the top and bottom of the wires. The ways to 
include these features into the thermal noise calcu- 
lations are different enough that they have led to 
some confusion among the gravitational wave com- 
munity. Also, attention has been paid mostly to 



the horizontal motion of the suspension, although 
all modes (angular, transverse, and vertical) appear 
to some degree into the detector's noise. We present 
a method to calculate thermal noise that allows the 
prediction of the suspension thermal noise in all its 
6 degrees of freedom, from the energy dissipation 
due to the elasticity of the suspension wires. We 
also show how the contributions of thermal noise in 
different directions can be sensed by the interferom- 
eter through the laser beam position and direction. 
The results will follow from the consideration of the 
coupled equations of the suspension and the contin- 
uous wire, first presented in || for just the horizon- 
tal degree of freedom. We show how this approach 
encompasses and explains previous ways to approx- 
imate the thermal noise limit in gravitational wave 
detectors. We show how this approach can be ex- 
tended to more complicated suspensions to be used 
in future LIGO detectors. To our knowledge, this is 
the first time the thermal noise of angular degrees of 
freedom is presented, and that all suspension degrees 
of freedom are calculated in an unified approach. 

Since the full treatment of the problem is some- 
what involved, we present first the problem without 
considering the elasticity of the wire, but adding a 
second, lossy, energy source to the gravitational en- 
ergy in the treatment of the mechanical pendulum, 
and introduce the concepts of "dilution factors" , and 
"effective" quality factors. We also start with one 
and two-degrees of freedom suspensions instead of 
6-dof. With these tools, most of the issues can be 
clearly presented and then we follow to the full treat- 
ment of the LIGO suspensions, presenting the impli- 
cations for LIGO. 
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II. SIMPLE PENDULUM CASES: DILUTION 
FACTORS, COUPLED MODES, EFFECTIVE 
QUALITY FACTORS. 



B. A pendulum with two energy sources: the 
dilution factor. 



The full treatment of this case, considering the 
elastic coupling of the wire to the suspension, was 
presented in [pj. Here, we will present the simpler 
"mechanical" treatment of this case, which will in- 
troduce the concepts of "dilution factors" , and mea- 
sured vs. effective quality factors. 

A. A simple oscillator with a dissipative energy 
source 



Next, we consider a suspended point mass, but 
we now assume there two sources of energy, gravi- 
tational and elastic, each with its own spring con- 
stant. The potential energy is then V = V g + V e = 
(l/2)(K g + K e )x 2 , and 

jmyl = u{K g (j> g + K e 4>e) 

' [ J ((K g + K e ) - Mlu 2 ) 2 + (K g + K,) 2 ^ 2 



We first recapitulate the calculation of thermal 
noise in the simplest case, a suspended point mass. 
The potential energy is PE = (l/2)Kx 2 and Fk — 
—dV/dx = —Kx. The kinetic energy is KE = 
(l/2)Mx 2 . The admittance to an external force F ext 
is given by 



Y 



K - Mlo 2 



The admittance has a pole at the system eigen- 
frequency w 2 , = K/M. If K is real, the resonance 
has an infinite amplitude and zero width. If the 
spring constant has an imaginary part representing 
an energy loss, K — > K(l + id>), then the amplitude 
|Y(w )| is finite, and the peak has a width deter- 
mined by the complex part of the eigenfrequency 
lo 2 = [K/M)(l + i<j)). The width of the peak is char- 
acterized with a quality factor Q — 1/0, and it is 
usually measured from the free decay time r of the 
natural oscillation at the frequency loq\ Q = loqt/2. 

The thermal noise is proportional to the real part 
of the admittance, and thus to dr. 



K[Y] = 



wKi> 



(K - Mu 2 ) 2 + K 2 (f> 2 

We are usually interested in frequencies well above 
too, since the pendulum frequency wo/2ir in gravita- 
tional wave detectors is usually below 1 Hz, and the 
detectors have their maximum sensitivity at 100 Hz. 
At those frequencies, the thermal noise is 

^BT 0mtx „ 4k B T Lo 2 cf> 



x 2 (uj) = 



-3t[Y] 



Mlo 5 



(1) 



This how we see that the measured decay of the 
pendulum mode can be used to predict the suspen- 
sion thermal noise at gravitational wave frequencies. 
Some beautiful examples of these difficult measure- 
ments and their use for gravitational wave detectors 
are presented in [§J, for example. 



If we assume that K g ^> K e , then 

LO(Kg(j) g + K e (j) e ) 



»ry] 



(K g - Mlo 2 ) 2 + K 2 ^ 2 
and at high frequencies 



x 2 (uj) - 4fcgJ o sft[y] ^ 4k B To(Kg<f> g + K e (j) e ) 



M 2 LJ 5 



If 



Kg(j> g < K, 



g = ("gravity is lossless"), or at least 
then 



where 



Ak B T uj 2 {K g /K e 
Mlo 5 



Kg/M. We see that is the same ex- 



pression as if we had just one energy source with a 
complex spring constant K = K g (l + i(K e / K g )d> e ). 
This is why we call the factor K e /K g the "dilution 
factor": the elastic energy is the one contributing 
the loss factor to the otherwise loss- free K g , but 
"diluted" by the small factor K e /K g . The dilution 
factor is also equal to the ratio of elastic energy to 
gravitational energy K e /K g = V e /V g . The concept 
of a dilution factor is very useful because it is usually 
easier to measure the loss factor d> e associated with 
the elastic spring constant than the quality factor of 
the pendulum mode. This is because K e is usually 
a function of the complex Young modulus E, and 
the imaginary part of the Young modulus is easily 
measurable for most fiber materials, and can even be 
found in tables of material properties. (Of course, 
there are subtleties to this argument, in particular 
with thermoelastic or surface losses ||, but we are 
assuming the minimum material loss). 
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C. A point mass suspended from an anelastic 
wire: calculating the dilution factor 

This case is a particular case of the one treated in 
||, and here we just mention it to present the ap- 
proach taken to the full problem, and present some 
new relevant aspects. 

We want to include the elasticity of the wire in 
the equations of motion, so we treat the suspension 
wire as an elastic beam, and then we have pendulum 
degree of freedom x, plus the wire's infinite degrees 
of freedom w(s) of transverse motion. We define a 
coordinate s, that starts at the top of the wire s = 0, 
and ends at the attachment point to the mirror, 
s = L. Correspondingly, we will have an eigenfre- 
quency, associated with the pendulum mode, and an 
infinite series of "violin" modes. The potential en- 
ergy is PE = (1/2) f p(w'(s)) 2 ds, and the kinetic 

energy is KE = (1/2) f£ p(w(s)) 2 ds + {l/2)Mx 2 . 
The solutions to the wire equation of motion, with 
boundary conditions w(0) = and w(L) = x 
are w(s) — xs'm(ks)/ sm(kL), and the equation 
of motion for the mass M subject to an external 
force F is F = -Mu 2 x + Tw'{L) = -Mw 2 x + 
(T/L)xQeL/taa(kL)). The admittance 

Y_. x _■ T tanfcL/fcL 

- IUJ— - ^J jl _ uj2 (ML Jjrj ( t£m kL / kL } 

has a pole at the pendulum frequency w 2 ~ T /ML, 
where kL <C 1, and an infinite number of poles at 
the violin mode frequencies, at frequencies w n ~ 
(T/pL 2 ) 1/2 (l + (pL/M)n 2 n 2 ), where kL = rnr + 
y/(pL/M)/(mr). 

The spring "constant" associated with the 
wire and gravity's restoring force K = 
(T/L)(kL/ tan(feX)) is in fact a function of fre- 
quency, although it is the usual constant T/L 
for frequencies below the violin modes, where 
kL <C 1. At frequencies above the first vio- 
lin mode, the spring function is not even posi- 
tive definite, or finite. The function K is real 
at all frequencies because we haven't added any 
source of energy loss yet. We introduce energy 
loss in the system by adding the wire elastic en- 
ergy to the system, and then assuming a com- 
plex Young modulus. The potential energy is now 
PE = (T/2)(J L w' 2 (s)ds) + (EI/2)J L w" 2 (s)ds. 
The equation of motion for the wire is 



a fourth order equation with boundary conditions 
w(0) = 0, w'(0) = 0, and w(L) = x. The wire slope 
at the bottom, w'(L), is a free parameter (since we 
are assuming a point mass), and the variation of 
the Lagrangian with respect to w'(L) provides the 
fourth boundary condition for the wire: w"(L) = 0. 
We can find an exact solution for the wire shape as a 
function of x, trigonometric functions of fcs, and hy- 
perbolic functions of k e s, where k, k e are solutions to 
Tk 2 — EIk 4 + pui 2 = which approximate at low fre- 
quencies the perfect string wavenumber, k 2 ~ pu 2 /T 
and a constant "elastic" wavenumber , k 2 ~ T/EI 
j|. The distance A = y/EI/T is the characteristic 
elastic distance over which the wire bends, especially 
at top and bottom clamps. In LIGO test mass sus- 
pensions, A ~2mm, a small fraction of L = 0.45m. 
The approximations k 2 ~ puj 2 /T and k 2 ~ T/EI 
are valid for frequencies that satisfy ui 2 <C T 2 /4EIp, 
about 12 kHz for LIGO, so we will use them in the 
remainder of this article. It is also equivalent to 
fcA < 1. 

We also use an approximate solution for the wire 
shape, good to order e~ L l A {~ 10~" (!) for LIGO): 

w{s) = A (sin(fcs) — fcA cos(fcs)) 

+kA(Ae- s/A + Be- {L - S)/A ) (2) 



The coefficients A, B are functions of x and fc, and 
thus, functions of frequency: 



A 



sin(fcL) — fcAcos(fcL) 
B = xkA. 

In the limit A — > 0, we recover the perfect 
wire solution, w(s) — xsm(ks)/sin(kL). The ra- 
tio B / A measures how much more (or less) the wire 
bends at the bottom than at the top. The elas- 
tic energy is well approximated by the contribu- 
tion of the exponential terms in the wire shape, at 
top and bottom: PE e = (1/2) J EI(w"(s)) 2 ds = 
{1/2)TA 2 f{w"{s)) 2 ds - (l/2)Tk 2 L{A 2 + B 2 ). At 
low frequencies where kL <C 1, the ratio B/A ~ 
A/L <C 1, indicating that the wire bends much more 
at the top than at the bottom (recall this is a point 
mass). 

The equation of motion for the mass when there 
is an external force F is 



F 



T 



d 2 w(s) 
ds 2 



EI 



d w(s) 
ds 4 



+ puj 2 w(s) = 0, 



-Muj 2 x + Tw'(L) 
-Muj 2 x + Tw'{L) 
( kT 



-Muj z x + 



EIw"'{L) 



1 + 2fcA tan kL 



\tanfcL 1 - (fcA)/(tanfcL) 



x (3) 
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The ratio of the elastic force to the gravitational 
force, EIw"'(L)/Tw'(L), is of order fcA < 1, and 
thus it was dropped. If we now consider A complex, 
then the spring function 

kT 1 + 2k A tan kL 
~ tan kL 1 - (kA)/ (tan kL) 

is also complex, and the admittance will have a non- 
zero real part. At frequencies below the violin modes 
where kL -C 1, we have K ~ (T/L)(l + A/L), an 
expression that suggests a split between a real grav- 
itational spring constant K g = T/L and a complex 
elastic spring constant K e = TA/L 2 . However, this 
distinction can only be done in the approximation 
A/L <C 1, and low frequencies kL <C 1. In general, 
however, we cannot strictly derive separate gravi- 
tational and elastic spring constants from their re- 
spective potential energy expressions: notice that 
the total, complex spring constant ~ T/(L — A) was 
derived from the variation of the gravitational poten- 
tial energy term, which becomes complex because we 
use a wire shape involving the complex distance A, 
satisfying the boundary conditions. 

Where the approximations A/L -C 1, kL -C 1 
are valid, we can consider the case of two separate 
spring constants and thus a "dilution factor" for the 
pendulum loss, K e /K g = A/L ~ 1/232, where the 
numerical value corresponds to LIGO parameters in 
Appendix 1. However, if we numerically calculate 
the exact pendulum mode quality factor, we get 
1/Qp = 461/(/>: this would be the measured Q from 
a decay time of the pendulum mode, if there are no 
extra losses. Did we make a "factor of 2" mistake? 
In fact, this factor of 2 has haunted some people in 
the community (including myself) JTof , but there is a 
simple explanation. Since the elastic complex spring 
constant is proportional to A, and A is proportional 
to the square root of the Young modulus, then when 
we make E complex E — * E(l + icf>), we get K e ~ 
TA/L 2 -> [TA/L 2 )(1 +i(j)/2). That is, we get an 
extra dilution factor of two between the wire loss <f> 
and the pendulum loss: 4> p = (K g / K e )^sK e /'3tK e i = 
(K g /K e i)<f>/2 = (A/2L)<j) ~ 0/464, very close to the 
actual value. This teaches us that if the spring con- 
stant of the dissipative force is not just proportional 
to E, we will get correction factors d\ogK e /d\ogE. 

The thermal noise below the violin modes is well 
approximated by the thermal noise of a simple os- 
cillator, as in Eqn.^, with natural eigenfrequency 
w a = 9/L an d loss 4>q = A</>/(2L). We then call 
4>o = A<f)/(2L) the "effective" loss, in this case equal 
to the pendulum loss (but we will see this is not 
always the case). 



We saw that the complex spring constant T/(A — 
L)was split into gravitational T/L and elastic 
TA/L 2 components. However both were derived 
from the gravitational force F g = —Tw'(L), since 
the elastic force, F e i — EIw"'(L), was negligible. 
It is because the wire shape w(s) is different due 
to elasticity, that the function Tw'(L) is different 
from the pure gravitational expression Tx/L. The 
way we split gravitational and elastic contributions 
to the spring constant and then got a dilution factor, 
is only valid at low frequencies. So the argument we 
posed in the previous section about a dilution factor 
applied to the calculation in the thermal noise in the 
gravitational wave band is in priciple not applicable 
here, especially when taking into account that the 
total force was contributed by the variation of just 
the gravitational potential energy, with the elastic- 
ity in the wire shape. However, using the wire shape 
without low frequency approximations, we can nu- 
merically evaluate the integrals that make up the po- 
tential and elastic energies (using a real A), and com- 
pare the ratio V e i/V g = fT(w') 2 ds/ J EI(w") 2 ds 
with the "low frequency" dilution factor A/L. We 
show the calculation of elastic and gravitational po- 
tential energies, and their ratio, in Fig.|l|. At low fre- 
quencies, the ratio is constant, and equal to A/2L: 
this is the dilution factor between the wire loss <f> 
and the pendulum loss <j) p , also the one to use for a 
simple-oscillator approximation of the thermal noise. 
It is not the ratio of the "gravitational" and "elas- 
tic" spring constants at low frequencies, but as we 
explained, we had no reason to expect that, since 
V e /V g ^ K e /K g . At higher frequencies, the ra- 
tio Vei/V g is not constant, and it gives correctly 
the dilution factors for the quality factors of the 
violin modes. Notice that the loss at the violin 
modes increases with mode number, as noted in 0]: 
</> n = (A/L)(l + n 2 n 2 A/L) = (A/ X)(l + P LAuo 2 jT), 
and this anharmonic behavior is well followed by the 
energy ratio. 

In summary, the concept of dilution factor is 
strictly true only when the total restoring force can 
be split into two forces, one lossless and one dissi- 
pative, both represented with spring constants. In 
the general case, if we can only split the potential 
energy into two terms, one lossless and another dis- 
sipative, then the ratio of the energies calculated as 
a function of driving frequency is the exact dilution 
"factor" . Moreover, this ratio can be calculated as a 
function of frequency, and then we get the different 
dilution factors for all the modes in the system. This 
is an important lesson that also we will use more ex- 
tensively in suspensions with more coupled degrees 
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of freedom. 



D. A suspended extended mass: coupled 
degrees of freedom and observed thermal noise 

We now consider an extended mass instead of a 
point mass, with a single generic dissipative energy 
source. The pendulum motion is described with the 
horizontal displacement of its center of mass x, and 
the pitch angle of the mass, 9, as in the side view 
of a LIGO test mass, in Figjj. The kinetic energy is 
KE = (l/2)(Mx 2 + J9 2 ). Instead of a spring con- 
stant, we have a 2x2 spring matrix. The potential 
energy is 

T , 



PE= ^(La 2 + h9 2 ) 

^(K xx x 2 + 2K x0 x9 



(4) 



where a = (x + hO) / L and 9 are the normal coordi- 
nates. The point x + h6 is the point where the wire is 
attached to the mirror, and the angle a is the angle 
the wire makes with the vertical. If we only consider 
gravitational forces, K xx = T/L, K x g = Th/L and 
Keg — Th(L + h)/L. However, we will assume that 
the elements of the spring matrix can be complex, 
and each has its own different imaginary part. The 
eigenfrequencies are the solutions to the equation 
(K xx - Muj 2 ){Kgg - Juj 2 ) - K 2 e = 0, or 



2MJw% = {MK ee + JK XX ) 
± ({MKee - JK X 



AMJKi 



.1/2 



(5) 



In order to calculate x 2 (to) (the Brownian motion 
of the center of mass), we need to calculate the ad- 
mittance x/F to a horizontal force applied at the 
center of mass. In order to calculate 9 2 (lo), we need 
the admittance 9/N to a torque applied around the 
pitch axis. If the spring constants are complex, then 
the admittances are complex and we can calculate 
their real parts, and the thermal noise determined 
by them: 



2 4k B To 

X (U!) = —M 



2 4fc B T 
1 (lj) = j — ^ 



1 - Juj 2 /K e 



K xx (1~co 2 /cuI)(1-u 2 /uj 2 _) 



1 - Muj 2 /K xx 

iMi-^ 2 M)(i-^ 2 /^) 



where the eigenfrequencies are now complex: ui\ — > 
+ i<t>±)- The quality factors measurable from 



the free decay of each of the eigenfrequencies are 
Q± = w±r±/2. 

At frequencies larger than any of the eigenfrequen- 
cies, we obtain 



2 ' - ^^((u 2 + ct> + +u J 2 _^)~c! > ee(Kee/J)) 
K xx uj° 

~ M((^ ++w 2^) _ ^ xx{Kxx /m)) (6) 



x 2 (lo) 



The system may have the two eigenfrequencies 
close in value if h(L + h) ~ J/M (see Fig.||), but 
for h < J /ML, we have cut — K ee /J = Th/J and 
w+ - K XX /ML 2 = T/ML; and for h(L+h) > J/M, 
uj\ - Kqq/J = Th/J and uj 2 _ - K XX /ML 2 = 
T/ML. In both limits, two terms cancel in the 
sum of loss factors in the formulas above {u) 2 _§- ~ 
(peeKee/ J for small h, for example) and we see that 



x 2 {u) 



4k B T uj 2 (f> xx 
Mw 5 



and 9 2 {u) 



4k B T uj^(t>e 
Juj 5 



Thus, even though it is a coupled system, ther- 
mal noise in x is always associated mostly to 
K xx (l + i<p xx ) and the pendulum eigenfrequency 
uo 2 = K xx /M ~ T/ML; and thermal noise in 9 is 
always associated with Kee(l + ifiee) and the pitch 
eigenfrequency uig — K 9 e/J ~ Th/J. For both de- 
grees of freedom x and 9, we obtain the thermal 
noise of single-dof systems. Unfortunately, neither 
limit (small or large h with respect to J /ML) applies 
to the suspension parameters in LIGO test masses, 
and, more importantly, even though the approxima- 
tion for the eigenfrequencies is relatively good for 
most values of h, the approximation we used for the 
losses is not (Fig.||) . The measurable quality factors 
give us <fi±, but we need 4> xx and 4>ee to use in the 
thermal noise of the pendulum, and these cannot be 
precisely calculated from <f>± unless we know 4> x e, or 
a way to relate it to the other loss factors. We will 
do this in the next section, using the elasticity of the 
wire. 

Notice that the forces and torques we have used 
to calculate the admittances Y x and Yg, will each 
produce both displacement and rotation of the pen- 
dulum. This means that the thermal noise in dis- 
placement and angle are not uncorrelated. This can 
be exploited to find a point other than the center 
of mass where the laser beam in the interferome- 
ter would be sensing less displacement thermal noise 
than at the center of the mirror, as was done follow- 
ing a somewhat different logic in ||. If we were 
to calculate the thermal noise at a point a distance 
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d above the center of mass, we then need to cal- 
culate the admittance of the velocity of that point 
(iuj(x + dO) = iu>x) to a horizontal force applied at 
that point. The equations of motion are 

F = (K xx - Mlo 2 )x + K xB 9 

Fd = K x gX + [Kgg - J y U 2 )9 

and then the thermal noise is 



x 2 M 



4fcs7o 

u 2 
4fc B T 

UJ 2 



5R 



(x + dff) 



F 

- d 2 Yg, 



2dY x , 



= x 2 (uj) + d 2 9 2 (uj) + 2d^^3t(Y xe ) (7) 

where Y xx is the admittance of x to a pure force F, 
Ygg is the admittance of 9 to a pure torque N, and 
Y x g is the admittance of a displacement s to a pure 
torque N, equal to the admittance of 9 to a pure 
force F. There is an optimal distance d below the 
center of mass for which the thermal noise x 2 ^) is 
a minimum: this distance is d = —3l(Y x g) /^R(Ygg) . 
The resulting thermal noise is 



XminM = x 2 (u)) 



4fc B T (jfty; e ) 2 

W 2 ®Y 99 



which is less than the thermal noise x 2 (uj) observed 
at the center of mass. However, the expression ob- 
tained for the distance d is frequency-dependent: 
that means we have to choose a frequency at which 
to optimize the sampling point. 

Summarizing, we have shown that whenever there 
are coupled motions, the thermal noise sensed at a 
point whose position depends on both coordinates is 
not the sum in quadrature of the two thermal noise 
[x 2 and d 2 9 2 {oj) in our case), but a combination that 
depends on the "cross-admittance". Moreover, the 
thermal noise of each degree of freedom cannot in 
general be calculated just from the measured qual- 
ity factors if the modes are coupled to each other 
strongly enough. 



E. A 2-DOF pendulum suspended from a 
continuum wire 

We add to the previous 2-DOF pendulum a con- 
tinuum wire, to be able to add the losses due to the 
wire's elasticity, and calculate modal and effective 
quality factors, as well as the point on the mirror at 
which wc can sense the minimum thermal noise. 



If we add elasticity to the problem, as in S, the 
potential energy is PE = {T/2)(J^ w' 2 (s)ds+h6 2 ) + 

(EI/2) w" 2 (s)ds. The boundary conditions for 
the wire equation are w(0) = 0, w(L) — x + hd, and 
w'(0) = 0, w'(L) = —9. The equations of motion 
for the pendulum are 



- Muj 2 x + Tw'(L) - EIw"'(L) = F 
-Joj 2 9 + EI{w"(L) + hw"'(L)) = N 



(8) 



In order to complete the equations of motion of 
the pendulum, we need the shape of the wire at the 
bottom end. For this, we use the shape given by 
the expression in Eqn.^, but this time the top and 
bottom weights are given by 



.4 



(h + A)9 



D 



B = ——(kx(cos(kL) 
kD 

+6»(sin(fc J L) 



fcAsin(fcL)) 
• k(h - A) cos(fci)) 



with D = sin(fcL) - 2kA cos(fcL) and A = y/EI/T 
as we used earlier. With the shape known, we can 
write the equations for x, 9 with a spring matrix: 



Ka T X 



Muj 2 )x + K x g9 
-(Kg 



JyLO 



F 

N 



where the spring functions are 



K xx = Tk(cos(kL) + kAsin(kL))/D 

Kgg = T(h + A)(sin(fcL) + k(h - A) cos(kL))/D 

K x g = K ex = Tk(h + A)(cos(Jbi) + kAsin(kL))/D (9) 

At this point, even though the expressions are 
complicated, we can calculate the complex admit- 
tances Y xx — ioJx/F, Ygg = iu9/N using a com- 
plex E and A. The analytical expressions for the 
admittances are quite involved, but we can always 
calculate numerically the thermal noise associated 
with any set of parameters. We can also calculate 
the widths of the peaks in the admittance, which 
would correspond to measurable quality factors for 
the pendulum, pitch, and violin modes. The plots 
presented in Figs, g and || were calculated using 
these solutions. 

Fig.^| shows that the frequency and quality factor 
of the pendulum and pitch modes vary significantly 
with the pitch distance h. At frequencies close to the 
pendulum eigenfrequencies, the thermal noise spec- 
tral densities show peaks at both frequencies. How- 
ever, at higher frequencies, the thermal noise x 2 (uj) 
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can always be approximated by the thermal noise 
of a simple oscilaltor as in Eqn.|l], with an "effec- 
tive" quality factor that fits the amplitude at high 
frequencies to the position of the single peak. We 
can similarly define an effective quality factor for 
the thermal noise in 9 2 (lo). We show in Fig.|| the 
actual and approximated thermal noise, with their 
corresponding effective quality factors found to fit 
best at 50 Hz. The effective quality factor at 50 
Hz can be calculated as a function of pitch distance, 
and we show this calculation in Fig. 0. The effec- 
tive pitch quality factor is well approximated, for 
any pitch distance, by the measurable pitch quality 
factor, while the pendulum effective quality factor is 
close to the measurable quality factor of the pendu- 
lum mode only at very small, or very large pitch 
distances. For the L1GO pitch distance of 8mm, 
the measurable pendulum quality factor is f times 
lower than the effective quality factor, and would 
then give a pessimistic estimate of thermal noise am- 
plitude. 

Low frequency approximation. At low frequencies, 
where kL <C f , we have expressions that can help us 
understand how the elasticity loss factor contributes 
to the effective quality factors, as well as to the pen- 
dulum and pitch modes. We trade this gain in sim- 
plicity for the loss of expressions valid at or above 
violin mode resonances. 

The low frequency limit of the spring constants in 
Eqns. ^| is 

K xx =T/(i-2A) 

K ee = T(h + A)(L + h- A)/(L - 2A) 
K 0X = K x6 = T(h + A)/(L - 2A). 

If we assume A has an imaginary part related 
to the material <j>: A — > A(I + i(j)/2), then we get 
complex spring constants. If we use these complex 
spring constants in Eq.|| we can calculate the loss 
factors of the pendulum and pitch mode, 4> p and 
(j>e; and if we use them in Eqns. |] we can get 
the effective quality factors. Using /i/I < 1 and 
A/L < f, we get Q xeS = 1/4> XX ~ A(f>/L and 
QecS = ^/4>ee ~ A<fr/2(h + A). This represents 
a "dilution factor" in displacement of A/L and in 
pitch of A/2(h + A). These approximations fit very 
well the values shown in Fig. 2|. 

At low frequencies, the equations of motion for 
x, 9 can be derived from a potential energy 



Using a complex A, this gives us the complex spring 



constants that we can use to get mode and effec- 
tive quality factors. We would like to break up 
this potential energy into a gravitational part and 
an elastic part, corresponding to a real gravitational 
spring constant (independent of A) and a lossy elas- 
tic spring constant. We know that if we take the 
limit A — > in P-Efci,<i, we obtain the regular po- 
tential energy PE g in Eq^ for a 2 DOF pendulum 
without elasticity. Thus, we are tempted to say that 
the elastic energy is the remainder, proportional to 
A, and thus having a complex spring constant when 
we consider a complex A. According to this argu- 
ment, we get K gxx = T/L and K exx = 2TA/L 2 . 
This would then give us a dilution factor for the dis- 
placement loss K e /K g ~ 2A/L. This is the factor 
by which the imaginary part of K e i is diluted; how- 
ever, as explained before, we pick up another factor 
of two due to K e being proportional to \f~E~. Thus, 
the dilution factor between the effective quality fac- 
tor and the wire quality factor is A/L. 

Energy Ratios and the Dilution Factor.We have 
seen that there is another way of identifying the 
"gravitational" and "elastic" terms in the poten- 
tial energy, using the actual potential energy expres- 
sions from which the equations of motion were de- 
rived: PE g = (T/2)(/ L w' 2 (s)ds + h8 2 ) and PE ei = 

(EI/2) w" 2 (s)ds. For any given applied force, or 
torque, we can solve the equations of motion for x, 9 
and w(s). (We don't need to invoke a low frequency 
approximation to do this calculation.) Then, if we 
calculate the ratio of elastic to gravitational poten- 
tial energy for a unit applied force, we get a func- 
tion which is frequency dependent, and is equal to 
the dilution factors for the pitch mode at the pitch 
frequency, for the pendulum mode at the pendulum 
frequency, and for the effective quality factor at high 
frequencies, A/L. We show the energy values and 
ratios for different values of the pitch distance in Fig. 

I 

Potential Energy Densities. There is another inter- 
esting calculation we can do with the solution ob- 
tained for the wire shape, and that is to find out 
where in the wire the elastic potential energy is con- 
centrated. In other words, we want to find a rela- 
tionship between the variation of the dilution fac- 
tor with frequency, and the curvature of the wire, 
mostly at the top and bottom clamps. Since both 
the gravitational energy and the elastic energy in- 
volve integrals over the wire length, we can define 
energy densities along the wire, and calculate a cu- 
mulative integral from top to bottom. The gravi- 
tational potential energy has also a term Th9 2 /2: 
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we define a ratio R = Th0 2 /(J w l2 ds), indicat- 
ing the relative contribution of this "pitch" term. 
From Fig.^, we observe that the gravitational po- 
tential energy density is distributed quite homoge- 
nously along the wire, even at the first violin mode. 
However, the pitch term, which can be considered 
a "bottom" contribution, contributes most of the 
gravitational energy when the system is excited at 
the pitch eigenfrequency, but also in several other 
cases at the pendulum frequency and at low frequen- 
cies. The elastic energy density is concentrated at 
top and bottom portions of length 2A. At low fre- 
quencies, the top contributes the most; at the pen- 
dulum eigenfrequency, the relative contributions de- 
pend strongly on h, but the bottom contributes at 
least half the energy; at the pitch eigenfrequency, the 
bottom contributes more than 99% of the energy; at 
higher frequencies, including the violin modes, top 
and bottom contribute equally. 

Motion of points away from center of mass. We 
discussed previously how it was possible to find a 
point whose thermal noise displacement was smaller 
than the thermal noise displacement of the center of 
mass. Now that we have expressions for the com- 
plex spring functions, we can find the optimal point 
and discuss the differences. The cross admittance in 
Eq.@ is 



_ K*e 

{K xx ~Mtu 2 )(K e9 - Ju 2 ) 

MJuj 4 

and then the optimal point (otimized at frequen- 
cies in the gravitational wave band, above pendulum 
modes) is d = -U(Y x e) /$t(Y ee ) ~ J /ML. Notice 
that even though the optimal distance was deduced 
from the thermal noise expressions, which all involve 
loss factors, the optimal distance only depends on 
mechanical parameters. As first explained in ||, the 
interpretation of this distance is that when the pen- 
dulum is pushed at that point by a horizontal force, 
the wire doesn't bend at the bottom clamp, produc- 
ing less losses. The fact that we can recover the 
result from the FDT is another manifestation of the 
deep relationship between thermal fluctuations and 
energy dissipation. We show in Fig^| the dependence 
of the thermal noise at 50 Hz on the point probed 
by the laser beam on the mirror, and the ratio of the 
thermal noise for d = and d = do at all frequen- 
cies. As expected, since the integrated rms has to be 
the the same for any distance at which we sense the 



motion, the fact that the spectral density is smaller 
at 50 Hz if d = do means that the noise will be 
increased at some other frequencies: this happens 
mainly at frequencies below the pendulum modes. 

There are many lessons to be learned from this ex- 
ercise, but perhaps the most important one is that 
the explicit solutions to the equations of motion have 
many different important results: 

• using the solutions to calculate the elastic and 
gravitational potential energies allows us to 
calculate a "dilution function" of frequency, 
equal to the dilution factor at each of the res- 
onant modes of the system, as well as to the 
most important effective dilution factor at fre- 
quencies in the gravitational wave band; 

• we can calculate energy densities along the 
wire to identify the portions of the wire most 
responsible for the energy loss and thus the 
thermal noise; 

• we can use low frequency approximations to 
find out expressions for dilution factors that 
can be found using other methods, explaining 
in this way subtleties like factors of two; 

• we can calculate the admittance of an arbi- 
trary point in the mirror surface to the driv- 
ing force, and thus find out improvements or 
degradation of observed noise due to beam 
misalignments. 



III. THE LIGO SUSPENSIONS: THERMAL 
NOISE OF ALL DEGREES OF FREEDOM 

We will now calculate the solutions to the equa- 
tions of motion for the six degrees of freedom of a 
LIGO suspended test mass, and then use the solu- 
tions to calculate the thermal nosie of all degrees of 
freedom, as well as the observed tehrmal noise in the 
gravitational wave detector. 

The mirrors at LIGO are suspended by a single 
wire looping around the cylindrical mass, attached 
at the top at a distance smaller than the mirror di- 
ameter, to provide a low yaw eigenfrequency. This 
is equivalent to having a mass suspended by two 
wires, attached slightly above the horizontal plane 
where the center of mass is. The mirror's 6 degrees 
of freedom are the longitudinal and transverse hor- 
izontal x and y, and the vertical z, displacements 
of the center of mass; the pitch 6 and yaw cj) ro- 
tations around the y and z axis, respectively, and 
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the roll ip around the longitudinal x axis. We show 
the coordinate system used and the relevant dimen- 
sions in Fig.^. The parameters used in the calcula- 
tion presented are those for LIGO test mass suspen- 
sions (Large Optics Suspensions). The mass of the 
cylindrical mirror is 10.3 Kg, the diameter is 25cm, 
and the thickness 10cm. The cylindrical wires are 
made of steel with density p = 7.8 x 10 3 kg/m 3 and 
0.62mm diameter. We assumed a complex Young 
modulus E = 2.1 X 10 1 1(1 + 10" 3 i)kg/m 2 . The ver- 
tical distance between the center of mass and the 
top clamps is I — 45cm, the wires are attached to 
the mass a distance h = 8.2mm above the center 
of mass. The distance between the top attachment 
points is 2a = 33.3mm. (In the previous examples 
where one wire was used, we assumed the same wire 
material and the same test mirror, but we used a 
0.88mm radius, so the stress in the wires remained 
constant.) 

Each wire element has displacement in a 2- 
dimensional plane transverse to the wire, Wi±(s) and 
a longitudinal displacement along the wire, w,ii(s). 
The kinetic energy is given by 

KE = \(T, f ' p\M»)?d* 

+M (x 2 + y 2 + z 2 ) + Jjj' 2 + J y 9 2 + J z 4> 2 ) 

The potential energy is given by the sum of the ax- 
ial strain energy and the bending (transverse) strain 
energy in each wire: 

PEt = lJ Q ds { Tw ^ 2 + EA ( w 'i\\f) 
+ l - J\sEI « ± ) 2 

plus the energy involved in rotating the mass: 

PE M = T {h cos a(9 2 +ip 2 ) - bsma{p 2 + ip 2 )). 

The wires will be attached at the top (s = 0) 
at the coordinates Wi — (0,±a, I), where / is the 
vertical distance of the top support from the equi- 
librium position of the center of mass. The wires' 
transverse slopes at the top will be zero. At the 
bottom the wires are clamped to the mass a dis- 
tance h above the center of mass, and a distance 2b 
on the y-direction between the wires on each side of 
the mass. The angle a = arctan((6 — a)/{l — h)) is 
the angle at which the wires are slanted from top to 
bottom when looking at the mass along the optical 
axis. If b — a, the wires hang vertically. The length 



of the wires is L 2 = (b — a) 2 + (I — h) 2 . The ten- 
sion in each wire is T = Mg cos a/2. The position 
of the bottom attachments when the mirror is mov- 
ing with a motion described by (x,y,z,6,<p,ip) are 
Wi (L) = (x + hO ± bp, y — hip, z±bip), and the slopes 
at the bottom are w[{L) = (—6, tp, 0). 

If we express the wire transverse and longitudi- 
nal displacements in the x,y,z coordinate system, 
we have Win = —w z cosa ± w^sina, and Wij_ = 

y/wi 2 , + (wi y cos a±Wi z sin a) 2 and the equations 
of motion become non-linear. In order to keep the 
problem simple, without losing any degree of free- 
dom, we will then consider two different cases: (i) 
the wire only has displacements in the x direction, 
and the mirror moves in x,9,4> degrees of freedom; 
and (ii) the wire only has displacements in the y, z 
directions, and the mirror moves in y, z, ip degrees of 
freedom. We analyze these cases separately. 

A. Longitudinal, Pitch and Yaw Thermal Noise 

The boundary conditions for the wires at the top 
are zero displacements and slopes, and at the bot- 
tom attachment to the mass, w±i(L) = x + h6 ± 

<j>,v/ ±i (L) = -e. 

We combine the wires' transverse displacements 
into w±(s) = (wi(s) ± W2{s))/2, then the bound- 
ary conditions at the bottom are w+(L) = x — 
h6,w' + {L) = 9,W-{L) = bp,w'_{L) = 0. The solu- 
tions to the wire equations of motion and the bound- 
ary conditions are (up to order e~ L ^ A ) will then be 
w±(s) as in Eqn.||, with 

A- = cf>b/D 

B_ = -<pb(cos(kL) + kA sm(kL))/D 

A+ = (x + (h + A)9)/D 

B + = -A(fcr(cos(fcL) + fcAsin(fcL)) 

+9(sm(kL) + k(h - A) cos(kL)))/D 

where D = sin(fcL) — 2fcAcos(fcL). The equations 
for the mass transverse dof subject to a force F x and 
torques N y ,N z , are 

Fx = -Muj 2 x - 2EIw'"{L) + 2T9 

N y = -J y u 2 9 + 2EI (w'l(L) + hw'"{L)) 

+2Th9(l - cos a) 
N z = -J z uo 2 4> - 2EIbw m {L) - 2Tb<j)sma (10) 

We see that the combination w + (s) is associated 
with the x, 9 degrees of freedom just as for the single 
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wire case, while the combination w- (s) is associated 
with the yaw degree of freedom <f>. This is easily un- 
derstood when imagining the wires moving back and 
forth "in phase" (io_ = 0), producing displacement 
and pitch but nt yaw; while if they move back and 
forth in opposition (w + = 0), then the only effect 
is into mirror's yaw. Thus, we can solve the equa- 
tions for w- (s) and <f> separately from the equations 
for uj + (s), x and 6: we will do so in the next parap- 
graphs. 

Yaw angular thermal noise. The admittance 
of yaw <f> to a torque N z is 

Ya, = iU ; 

— 2Tb sin a — J z uj 2 

with (2Tb 2 /L)(kL)(cos(kL) + 

kAsin(kL))/ (sin(kL) — 2fcAcos(fcL)). As usual, 
when A = y/EI/T is complex, the admittance is 
complex and using the fluctuation-dissipation theo- 
rem we obtain cf) 2 {f) = (4fc B T o )5R(y )/w 2 . 

The effect of the tilted wires with a — (b — a)/L > 
is to lower the restoring force, and thus the res- 
onance frequency: in LIGO suspensions, the fre- 
quency is 0.48 Hz instead of 1.32 Hz if b = a. 
However, since sin a ^ decreases the gravita- 
tional restoring force but not the elastic force, the 
dilution factor ~ K e /K g increases, and so docs 
the thermal noise. At frequencies where kL <C 1, 
~ (2Tb 2 /L)(l + 2A/L). The thermal noise at 
frequencies below the violin modes, where kL <C 1, 
is well approximated by the thermal noise of a single 
oscillator with resonance frequency = 2Tab/JL 
and quality factor Q$ = (L/A)(a/b)(l/(f> w ), where 
E — ► E(l + icf> w ). Thus, the "dilution factor" is 
(A/L)(b/a) (=1/72 for LIGO parameters, where 
b/a = 7.5). As in the case of a simple pendulum, 
this dilution factor is half of the ratio of the elastic 
spring constant K e = (2Tb 2 / L)(2A/ L) to the grav- 
itational spring constant K g = 2Tab/ L, because the 
elastic spring constant has an "extra" dilution factor 
of 2: if E -> E(l + i<f> w ), then A — > A(l + i(j) w /2) 
and K e — > K e (l+i<j> w /2). However, the ratio of elas- 
tic potential energy J L EI(w±) 2 ds to gravitational 

potential energy V g — T(w' ± ) 2 ds — Tbsinacj) 2 is 
equal to the "right" dilution factor (b/a) (AIL) ~ 
1/44 at low frequencies, as shown in Fig. |[ The 
energy ratio also gives us the right dilution factor at 
the violin frequencies (2A/L). 

The yaw angular thermal noise may be seen in 
the detectors' gravitational wave signal if the beam 
hits a mirror at distance d to either side of the cen- 
ter of mass, or if it hits the mirror in a direction 



an angle 7 away from longitudinal. Considering 
both cases, the sensed thermal noise will be given 
by x 2 (w) = (dcos7 + H sin7/2) 2 2 (w), where H is 
the thickness of the mirror: 

X 2 H = (dcos 7 + (H/2)sm 1 ) 2 <j) 2 (Lj) 

AksT , , TTI . . , 2 2T6 2 A , 

\- (d cos 7 + (H/2) sin 7 ) 2 — - - 

J X L L 

where the approximation is valid between the pendu- 
lum mode and the first violin mode, lHz-50Hz. At 
160 Hz, where the maximum sensitivity of \fh 2 (f) = 
2.5 x 10~ 23 / ^VHz is expected, the yaw thermal noise 
is y<F(160Hz) = 2.9 x 10~ 19 rad/\/Hi. If the yaw 
thermal nose is to be kept an order of magnitude 
below the dominant noise source, then it is required 
that d <lcm and 7 < 10°. 

Notice that the mirror will always be aligned nor- 
mal to the laser beam to make the optical cavities 
resonant; however, what matters is the beam direc- 
tion with respect to the coordinate system defined by 
the local vertical and the plane defined by the mirror 
in equilibrium. Presumably there will be forces ap- 
plied to align the mirror, but in principle they have 
no effect on the response of the mirror to an osculla- 
tory driving force such as the one we imagine in the 
beam's direction, to calculate the admittance. Thus 
the requirement on 7 < 10° is on the position of the 
mirror when there are no bias forces acting, with re- 
spect to the ultimate direction of the beam. The 
beam's direction must be within 1/irad of the nor- 
mal to the aligned mirror to keep the beam aligned 
on mirrors 4km apart, but that doesn't mean that 
the mirror has not been biased by less than 10° to 
get it to the final position. 

Pitch and displacement thermal noise. We 
now solve the equations for w+(s),9 and x. The 
equations for these degrees of freedom in Eqns. [lC] 
are exactly the same as for the pendulum suspended 
on a single wire (Eqn. ||), except for the addition of 
a softening term to the torque equation, due to the 
tilted wires; and factors of two due to the two wires 
(with about half the tension) instead of a single wire. 
The extra term in the torque equation is a negligible 
contribution to the real part of Kgg, at the level of 
1% for LIGO parameters. Therefore, the conclusions 
we obtained, with respect to the optimization of the 
beam location on the mirror, and the difference be- 
tween effective and measurable quality factors, are 
equally valid here. The spring constants we obtained 
in Eqns. || involve now a factor 2T — Mg cos a, in- 
stead of T = Mg for a single wire. The elastic dis- 
tance is however determined by the tension in each 
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wire, A = y/EI/T = y/2EI/Mgcosa. To keep the 
stress in the wires constant, the cross section area 
of a single supporting wire is twice the area of two 
each of two supporting wires. Thus, the effective 
quality factor determining the thermal noise for the 
displacement thermal noise x 2 (f) has a smaller di- 
lution factor of A/L = 1/326, instead of 1/231 for a 
single wire. This is the well-know effect of reducing 
thermal noise by increasing the number the wires. 
The dilution factor for pitch is 1/14, considerably 
higher than the dilution factor for yaw, 1/44. 

The displa cement thermal noise at 160 Hz is 
■ s /x 2 (160Bz) = 1.1 x 10" 2O m/VFh7, limiting the de- 
tector sensitivity to h = 5.6 x 10~ 24 /vTIz- This is 
expected to be lower than the thermal noise due to 
the internal modes of the mirror mass, not consid- 
ered here p"l{ ]. The pendulum thermal noise could 
be reduced by a factor v2j or about 40%, if the 
beam spot was positioned at the optimal position 
on the mirror. Since pendulum thermal noise is not 
the dominant source noise, but the detectors' shot 
noise would increase due to diffraction losses, it is 
not advisable for LIGO to proceed this way. How- 
ever, these considerations should be taken into ac- 
count for future detectors, where thermal noise may 
be a severe limitation at low frequencies. 

The pitch angular noise at 160 Hz, is \J 2 (f) = 
8.9 x 10~ 19 rad\/Hz. Its contribution to the sensed 
motion has to take into account the coupling with 
displacement, and we will do this in detail inthe last 
section. 



B. Vertical, transverse displacement and roll 

We are now concerned with the mirror 
motion in its y, z and ip degrees of free- 
dom. The potential energy is PEi = 

(1/2) J Q L (T(w' ± y + EI(w'> ± y + EA(w> ll )*) ds for 

each wire, plus T(ft,cosa — b sin a)ip 2 . Notice that 
due to the tilting of the wires, the "transverse" w± 
and "axial" w\\ directions are not y and z, but rota- 
tions of these directions by the wire tilt angle a. We 
define Wij_ for each wire pointing "out" (and thus 
in opposite directions if a = 0), and Win pointing 
down along the wire, from top to bottom. 

The boundary conditions at top are Wi±(0) = 
0,^ x (0) = 0, lOiii(O) = 0, and at the bottom, 
Wij_(L) = ±(ycosa — dip) + zsina, w' i± (L) = ±tp, 
and ittj|i(L) = ±(ysina — cip) — zcosa, where we 
defined two new distances c = h sin a + b cos a and 
d = hcosa ~ b sin a. If we define as earlier, sums 



and differences of the two wires shape functions, 
w±(s) = (wi(s) ± W2(s))/2, then the equations of 
motion for the mirror degrees of freedom, when sub- 
ject to external forces F y ,F z and a torque N x , are 

F y = -Muo 2 y + 2(Tw'_ x (L) - EIw'" ± (L)) cos a 
+2EAw'_., sin a 

F z = -Muj 2 z + 2(Tw' +± (L) - EIw^ ± (L))siaa 
— 2EAw' + ,, cos a 

N x = ~J x uj 2 ijj - 2EI(dw'" ± {L) + w'LJL)) 
+2EAcw'_ l{ {L) 

The solution for the wires' transverse motion 
w±j_(s) satisfying the boundary conditions up to or- 
der e _i / A and fcA are of the same form as in Eqn. 
^, with top and bottom weights equal to 

A_ = (y cos a-(d+ A)V0 /D 

B- = — ((y cosa — dip)(cos(kL) + fcAsin(fcL)) 

-ip($in(kL) - kAcos(kL))/k)/D 
A + = z sina/D 

B + = -z sin a(cos(kL) + fcAsin(fcL))/L» (11) 



The axial wire motion is 



i»±\\ = W ±\\( L ) 



sin(fc z s) 
' sin(fe z L) 



(12) 



where ui-^(L) = j/sina — cip and w+^(L) = 

—z cos a. The wavenumber functions are k 2 — 
puJ 2 /T, k\ = P lu 2 /EA. 

Even though the tilting of the wires produces more 
complicated formulas than in the pendulum-pitch- 
yaw case, the equations for the vertical motion de- 
couple from the equations from the transverse pen- 
dulum displacement and roll, similar to yaw decou- 
pling from pendulum and pitch. As before, if the 
wires move in phase (transverse or axially or both), 
they produce only vertical motion; but if they move 
in opposition, they produce side to side motion plus 
rotation around the optical axis. We analyze the two 
decoupled systems separately. 

Vertical thermal noise. Once we have solved 
the wire shape (from Eqnsp|JTl|, and |l2|), we can write 
the equation of motion of the wire vertical displace- 
ment as 



F 7 



-Muo 2 z+2 



(EA k z L 



\ L tan(fc z L) 



cos 2 a 



-Tk 



cos(kL) 



kAsin(kL) . 2 \ 

sin a z 



sin(fcL) - 2fcAcos(fcL) 
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or 

F z = -Mlu' 2 z + (K T sin 2 a + K E cos 2 a)z 

where we defined Kt as the spring constant that was 
used in the pendulum-pitch case, and K E ~ 2EA/L 
the spring constant of the wire. For LIGO parame- 
ters, and for usual wires, Kt/Ke <C 1. 

If the wires are not tilted and a = 0, we recover 
the simple case of vertical modes of a mirror hang- 
ing on a single wire. The restoring force is elastic 
and proportional to E, so there is no dilution fac- 
tor. The term added because of the wire tilting is 
a gravitational restoring force, much smaller than 
the elastic restoring force. Since it is also mostly 
real when considering a complex Young modulus E, 
it will not change significatively the loss terms, and 
thus the thermal noise. The wire tilting does add, 
however, the violin modes to the vertical motion, 
and it slightly decreases (by a factor cos a = 0.97) 
the frequency of the lowest vertical mode. 

The vertical thermal noise at 160 Hz is *J z 2 (/) = 
3.1 x 10~ 18 m/\/Hz, 260 times the pendulum ther- 
mal noise. This is due to the lower quality factor, 
and the higher mode frequency. However, vertical 
noise is sensed in the gravitational wave interferom- 
eter through the angle of the laser beam and the 
normal to the mirror surface, which is not less than 
the Earth's curvature over 4km, (0.6 mrad). At the 
minimum coupling (0.3 mrad for each mirror in the 
4km cavity), the contribution due to vertical ther- 
mal noise is 10% of the pendulum thermal noise. In 
advanced detectors, vertical modes are going to be at 
lower frequencies due to soft vertical supports, like 
in the suspensions used in the GEO600 interferom- 
eter, but the ratio of quality factors is just the me- 
chanical dilution factor, so the contribution of verti- 
cal thermal noise to sensed motion will be of order 
3xl0~ 4 x y/L/Axf z /f x ~ ZxlQ- i x(EL 2 /Mg) 1 ^, 
not necessarily a small number! 

Side pendulum and roll. The side motion and 
roll of the pendulum are not expected to appear 
in the interferometer signal, but it is usually the 
case that at least the high quality-factor resonances 
do appear through imperfect optic alignment. The 
equations for the system y, ip can be written as 

(Kyy - MU 2 )y - Ky^l/j = F y 

(Kyy - Jx^ 2 )^ - KipyV = N x 

with 

Kyy = Kt cos 2 a + Ke sin 2 a 
= K E c 2 + K T d(d + A) 



sm(kL) — fcA cos(kL) 
+ ^ + Km(kL)-2kAcos(kL) 
K y ^ = K^y = Kx{d + A) cos ol + Kec sin a (13) 

Within the (very good) approximation Kt/Ke <C 
1, we can prove that the eigenfrequencies are w 2 — 
Kt /M, fy = 0.75 Hz and lo 2 = K E b 2 /JJ 4 , = 19 
Hz. The corresponding loss factors are <p y ~ A<fi/2L 
and 4>^, ~ 4>. 

If the wires are perfectly straight, the thermal 
noise of the side-to-side pendulum motion below the 
violin modes Uo(f) is well approximated by that of 
a simple oscillator with eigenfrequency uj y and a di- 
lution factor A/2L. The thermal noise of the roll 
angular motion ip 2 (f) does not depend much on 
the wire tilt, and is well approximated (below vi- 
olin modes) by the thermal noise of a simple oscilla- 
tor with eigenfrequency w 2 = 2EA/J X L and quality 
factor equal to the free wire's quality factor (there 
is no dilution factor). 

For any small wire tilt, however, the spring con- 
stant Kyy in Eqn.^ has a large contribution of 
Ke, and thus the fluctuations increase as y 2 (f) ~ 
Vp{f) cos 2 a + {J x /Mb) 2 ^ 2 {f) sin 2 a, as seen in Fig. 
0. Since the roll eigenfrequency is higher than the 
side pendulum eigenfrequency, and its quality factor 
is lower, this is a significant increase in the thermal 
noise spectral density tp 2 (f), about a factor of 100 
for LIGO parameters. However, the gravitationalw 
ave detectors are mostly immune to the roll degree 
of freedom, as we will see later. 

Also, the tilting of the wires introduces vio- 
lin modes harmonics of fn^ — \J\EAf p)/{2L) ~ 

7i x 6kHz apart from the usual harmonics fn = 
y/(T/p)/(2L) ~ n x 300Hz. 

The violin modes that are most visible in the roll 
thermal noise are the harmonics of y /f {EA/p)/(2L) 
(and are strictly the only ones present if the wires 
are not tilted). 



C. Violin Modes 

We have explored the relationship between quality 
factors of pendulum modes and the "effective" qual- 
ity factor needed to predict the thermal noise of any 
given degree of freedom (seen only in sensitive in- 
terferometers) . For the most important longitudinal 
motion, we've seen that the effective quality factor is 
approximately equal to Q W L/A, where Q w — 1/4> W 
is the quality factor of the free wire, related to the 
imaginary part of the Young modulus. Under ideal 
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conditions where the pendulum mode is far away 
from the pitch mode in frequency, its quality factor 
is close to the effective quality factor, but as we have 
seen, the errors may be as large as 50%. 

The violin modes, approximately equal to /„ = 
riy/T I ' p/2L, appear in the horizontal motion of the 
pendulum in both directions, along the optical axis 
and transverse to it. The violin modes show some 
anharmonicity, as pointed in 0],with the frequen- 
cies slightly higher than nyT/p/2L, and the quality 
factors degrading with mode number. The complex 
eigenfrequencies are the solutions to the equation 

kL = nn + arctan(2fcA/(l - (fcA) 2 )) 

with k 2 — pu> 2 /T. This means that if we measure 
the quality factors of violin modes, and they follow 
the predicted anharmonic behavior in both direc- 
tions, we can assume the losses are only limited by 
wire losses. We can then predict the thermal noise in 
the gravitational wave band with more confidence, 
having also consistency checks with the quality fac- 
tors measured at the pendulum modes. The thermal 
amplitude of the peaks at the violin modes follows a 
simple l// 5 / 2 law, corrected by the change in longi- 
tudinal Qs. We show all these features in Fig.p"2[ 

D. Total Pendulum Thermal Noise in LIGO I 

The right way to calculate the total thermal noise 
observed in the interferometer signal is to calcu- 
late the pendulum response to an applied oscillating 
force in the direction of the laser beam. The pen- 
dulum responds to a force in all its six degrees of 
freedom, but the motion we are sensitive to is the 
motion projected on the laser beam's direction. 

If the laser beam is horizontal, and its direction 
passes through the mirror center of mass, it will be 
sensitive to only longitudinal displacement and pitch 
motion. If the beam is not horizontal, and for ex- 
ample is tilted up or down by an angle 7, (but still 
going through the center of mass) , we imagine an ap- 
plied force F applied in the beam's direction, with 
a horizontal component F cos 7 and a vertical com- 
ponent _Fsin7. The motion we are interested in is 
x cos 7 +z sin 7, since it is the direction sensed by the 
laser beam. The admittance we need to calculate is 
then Y — iuj(x cos j + z sin / F = iuj((x/F x ) COS7 + 
(z/F z )shi7) = iw(Y x cos 2 7 + Y z sin 2 7), where Y x 
is the response of the horizontal displacement to a 
horizontal force and Y z is the vertical response to 
a vertical force. A force applied in the beam's di- 
rection, with magnitude Fq, will have components 



in all 3 axes, and torques around all 3 axes too: 
F = F a (a x ,a y ,a z ) and N = Fo(D^,Dq, £)</,). The 
motion we are interested is, in general, \ = a x x + 
a y y + a z z + Dg9 + D^<j> -I- D^ip. If the motion in 
all 6-DOF was uncoupled, then each dof responds to 
just one component of the force or the torque, and 
the admittance we need would be just the sum of the 
admittances, each weighted by the square of a fac- 
tor a,; or a distance Di. However, only z and </> are 
decoupled dof from the rest, and x, 9 and y, ip form 
two coupled systems, for which we need to solve the 
response to a forces and torques together. 

We define an admittance Y x g = iu>(x/N y ) as the 
admittance of displacement x to an applied torque 
N x , and Yg x = iui(9/F x ) as the admittance of pitch 
to an applied horizontal force , and similar quanti- 
ties Y y ^,Y^ t y. Then, the total thermal noise sensed 
by the laser beam is 

X 2 (f) = ^^^{a 2 x Y xx + a x Dg(Y x g + Yg x ) + D 2 Ygg 

+a 2 Y yy + a y D^(Y y ^ + Y^y) + D^Yjp^ 
+a 2 z Y zz + D 2 Y^) 

or 

X 2 (/) = a 2 x x 2 (f) + D 2 6 2 (f) + ax Dg^^^(Y x g + Yg x ) 

UJ 

+a 2 y y 2 (f)+D 2 i/j 2 (f) 

+a y D^, SftftW + Y-^y) 

+a 2 z z 2 (f)+D 2 <j> 2 (f) 

so it is a weighted sum in quadrature of the ther- 
mal noise of different degrees of freedom, plus some 
cross-terms. These terms may be negative, so it is 
possible to choose an optimal set of parameters to 
minimize the sensed motion, as shown in @]. The 
weighting factors and distances are 

a x = cos(jy) cos(j z ) 
a y = sin(7 a )cos(7 z ) 
a z = sin (7 2 ) 

Dg = Rsin(j z ) + dzcos(j y ) cos(7,z) 
£>0 = cos(7 z )(d y cos(jy) - i?sin(7j / )) 
A/, = d y sin(7 z ) + d z sin^) cos(7 z ) 

Typical distances d y , d z are few mm at most, and 
typical angles 7 y are in the order of microradians, 
since an angle of such magnitude produces displace- 
ments of the order of millimeters at the beam at the 
other end of the arm, 4km away. However, the angle 
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7z cannot be less than half the arm length divided by 
the curvature of Earth, or j z > 3 x 1CT 4 . We show 
in Figjl3| the thermal noise sensed by a beam with 
dy = dz — 5mm , -f z — 3 x 10~ 4 and j y = 1/irad. 



IV. RESULTS AND CONCLUSIONS 

We have shown several general results, and then 
calculated predicted thermal motions for LIGO sus- 
pensions. First, we showed that if there is a sin- 
gle oscillator with two sources of potential energy, 
one with a spring constant K g and a dominant real 
part, and the other with a complex spring constant 
K e , with a dominant loss factor, then the dilution 
factor K e /K g gives us the ratio between the oscil- 
lator's quality factor (determining its thermal noise 
spectral density) and the loss factor of K e . However, 
the elastic loss 4> e — ^sK e /^RK e might itself have a 
small dilution factor with respect to the loss factor 
of the Young modulus, for example if K e ~ yE, 
then there is a dilution factor of 1/2. 

We also show that when two or more degrees of 
freedom are coupled, the measurable quality factors 
at resonance may not be as useful to predict the 
"effective" quality factor used in the thermal noise 
spectral density, unless the eigenfrequencies are far 
from each other. 

We showed that by using approximations of or- 
der e~ L / A , we can easily obtain wire shapes and 
equations of motion for the 6 degrees of freedom 
of the mirror, as a function of applied oscillating 
forces and torques. We think that this method will 
be most useful when applied to multiple pendulum 
systems such as those used in GEO600 and planned 
for advanced LIGO detectors 0,0. However, even 
in the simple pendulum case this approach allows 
us to calculate the gravitational and elastic poten- 
tial energy as a linear energy density along the wire, 
and the total energy. We showed that, for an ap- 
plied horizontal force, the gravitational potential en- 
ergy ((1/2) fTw' 2 ds) is homogeneously distributed 
along the wires, while the elastic potential energy 
((1/2) J EI(w") 2 ds) is concentrated at the top and 
bottom, but in different proportions depending on 
the frequency of the applied force (Fig. [|. We also 
calculate the ratio of total elastic potential energy to 
gravitational energy using the solutions for the wire 
shape, and show that this function of frequency cor- 
responds to the dilution factor for the eigenmode 
loss factors as well as for the effective quality fac- 
tor that allows us to calculate the thermal noise at 
gravitational wave frequencies (see Figs.[l],[I|,||). 



Applying our calculation of wire shapes and equa- 
tions of motion that include elasticity to LIGO sus- 
pensions, we show in Figs [l^ and [Tl] the resulting 
spectral densities of displacement and angular de- 
grees of freedom of the mirror. More importantly, 
we show in FigfL3| the resulting contribution of pen- 
dulum thermal noise to the LIGO sensitivity curve, 
assuming small misalignments in the sensing laser 
beam (5mm away from center of mass, 1/irad away 
from horizontal). As expected, the displacement de- 
gree of freedom is the one that dominates the con- 
tribution, but pitch noise contributes significantly 
(29% at 100 Hz) if the beam is 5mm above center. 
But this is not added in quadrature to the horizontal 
noise, which makes 81%: the coupled displacement- 
pitch thermal motion makes up 99% of the total 
thermal noise. A misplacement below the center of 
mass will reduce the observed thermal noise, as first 
noted in Q . The contribution of yaw thermal noise 
(11% at 100Hz) is smaller than that of pitch, but 
very comparable. The contribution of vertical noise 
due to the 4km length of the interferometer is 8% 
at 100 Hz. The side and roll motions are coupled, 
but the roll contribution dominates (due to the large 
angle j z and the assumed d z = 5mm) and is 0.7%, 
much smaller than the contributions of pitch, yaw 
and vertical degrees of freedom. When added in 
quadrature, the total thermal noise is 23% higher 
than the contribution of just the horizontal thermal 
noise. However, if the vertical misplacement of the 
beam is 5mm below the center of mass, instead of 
above, the total contribution of thermal noise is 89% 
of the horizontal thermal noise of the center mass. 
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FIG. 1. Top figure: Gravitational and potential energies for a suspended point mass as a function of frequency 
/, when excited by a sinusoidal force wiht frequency /. Bottom figure: ratio of energies, considered as a "dilution 
factor". The stars indicate the factor (f)/Qi, where <j> is the wire loss factor and Qi is the eigenmode quality factor. 
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FIG. 2. Eigenfrequencies and associated loss factors for a pendulum suspended on one wire. Quality factors and 
corresponding eigenfrequencies are represented with the same symbol (points or stars). The lines in the top figure 
are given by expressions for the pendulum frequency f p = U g/l/(2ir) (solid) and pitch frequency f$ = y/Th/J / (27r) 
(dashed). The lines in the bottom figure correspond to the effective quality factors for displacement (solid) and for 
pitch (dashed). Approximate expressions for the effective quality factors are Q v ~ L/A<f> and Qg ~ 2(h + A)/A<f>, 
respectively. 
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FIG. 3. Displacement (top) and pitch (bottom) thermal noise of a single-wire pendulum, calculated without ap- 
proximations (solid line), and approximated with a single mode with an effective Q (dotted line). The pitch distance 
used was 8.3 mm 
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FIG. 4. At the top, we plot the elastic and gravitational potential energies for a pitch distance h=8.2mm. In the 
bottom figure, we plot the ratio of elastic and gravitational energy for three different values of the pitch distance. 
The stars represent the dilution factor of the pendulum mode, at the pendulum frequency; the circles are the dilution 
factors of the pitch mode, at the pitch eigenfrequency; and the solid horizontal line is the dilution factor of the 
effective quality factor at 50 Hz. 
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FIG. 5. The figure shows cumulative potential energies along the wire (top is s = 0, bottom is s = 0.45m), when 
the pendulum is excited by a force a the center of mass of frequency /. The energies are shown for three different 
pitch distances in each graph: h=2mm (dashed), h=8mm (thick points) and h=20mm (dotted). The energies are also 
plotted for five different frequencies: / = 0.3 Hz, below all eigenmodes; / = f p orpendulum frequency, approximately 
0.75 Hz; / = f g = y/Th/J / '(2n) = (0.46, 0.69, 1.12)Hz for h = (2, 8, 20)mm respectively; f=50Hz, in the gravitational 
wave band; and / = 319Hz, the first violin mode frequency. 
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FIG. 6. In the top figure, we plot the thermal noise amplitude spectral density (ASD) at 50 Hz, as a function of 
the distance d at which the laser beam is positioned. In the middle figure, we plot the thermal noise ASD times f 2 ^ 5 
obtained for d = and for the optimal distance for / = 50Hz. In the bottom figure, we plot the ratio of the two 
curves in the middle figure. 
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Front view Side view To P view 

FIG. 7. LIGO pendulum suspension, with a single wire loop, attached slightly above the center of mass and angled 
towards the center. 
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FIG. 8. Potential energies and dilution factor for LIGO yaw degree of freedom 
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FIG. 9. Thermal noise of pendulum side and roll motion, for different values of the distance between the wires 
at the top. The wires are tilted with respect to the vertical direction at an angle sin a — (b — a)/L, where L is the 
length of the wires. 
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FIG. 10. Thermal noise of the three mirror's displacement degrees of freedom. 
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FIG. 11. Thermal noise of the three mirror's angular degrees of freedom. 
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FIG. 12. In the top figure, we show the frequencies of the resonant violin modes, scaled by f n — nW T/rho/2L, 
for longitudinal and transverse modes. In the middle figure, we show the corresponding quality factors. The solid 
line is Qeff/2 ~ L/2A. In the bottom figure, we show the amplitude of longitudinal and transverse thermal noise 
at the violin modes. The solid line shows a l// 5 ' 2 fall off from the first peak. 
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